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A technique is described for the design of fuel-conservative guidance systems and is applied to a system that 
was flight tested on board NASA’s augmentor wing jet STOL research aircraft. An important operational 
feature of the system is its ability to rapidly synthesize fuel-efficient trajectories for a large set of initial aircraft 
positions, altitudes, and headings. This feature allows the aircraft to be flown efficiently under conditions of 
changing winds and air traffic control vectors. Rapid synthesis of fuel-efficient trajectories Is accomplished in 
the airborne computer by fast-time trajectory integration using a simplified dynamic performance model of the 
aircraft. This technique also ensures optimum flap deployment and, for powered-lift STOL aircraft, optimum 
transition to low-speed flight. Also included in the design is accurate prediction of touchdown time for use In 
four-dimensional guidance applications. Flight test results have demonstrated that the automatically synthesized 
trajectories produce significant fuel savings relative to manually flown conventional approaches. 
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Nomenclature 

= drag force, lb 

= distance of backward and forward 
integration, respectively 
= cruise distance, ft 

= length of ground track from initial to final 
position of aircraft, ft 
= energy, ft 
= energy rate, ft/s 

= maximum and minimum available energy 
rate, respectively, ft/s 

= perturbation state and control distribution 
matrices, respectively 
= acceleration of gravity, ft/s^ 

= final and initial ground headings of 
aircraft, deg 
= altitude, ft 

= final and initial altitudes of aircraft, 
respectively, ft 
= feedback gain matrix 
= lateral error and error rate feedback gains 
= lift force, lb 

= speed along ground track, ft/s 
= thrust force, lb 
= time, s 

= perturbation control vector 
= two-dimensional unit vectors 
= aircraft control vector 
= airspeed, ft/s or knots 
= final and initial airspeeds of aircraft, 
respectively, ft/s or knots 
= wind speed in direction of ground track, 
knots 

= aircraft weight, lb 
= perturbation state vector 
= final and initial X coordinates of aircraft, 
respectively, ft 

= final and initial Y coordinates of aircraft, 
respectively, ft 
= angle of attack, deg 
= ineritial flight-path angle, deg 
= aerodynamic flight-path angle, rad or deg 
= airspeed rate correction due to wind shear, 
ft/s 2 
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= crosstrack error, ft 
= crosstrack error rate, ft/s 
= flap angle, deg 
= maximum flap angle, deg 
= fraction of energy rale used for changing 
speed 

= command pitch angle, deg 
= vectored thrust angle, in degrees of nozzle 
angle 

= throttle setting, in percent rpm 
= fraction of available energy rate 
= commanded and reference bank angles, 
respectively, deg 

Introduction 

I N the past, terminal area guidance system design for 
aircraft has concentrated primarily on automatic glide- 
slope tracking, flare, and touchdown. During recent years, 
designs have been developed to provide automatic guidance 
along curved and decelerating approach paths. ^ This in- 
creased capability was made possible through the integration 
of digital computers into the flight guidance system. 
However, even in the more advanced designs, automatic 
guidance is limited to a few prestored three-dimensional (3-D) 
flight paths, as in Ref. 1. While the ability to fly complex 
prestored trajectories is essential, it cannot give optimum 
performance under actual terminal area operating conditions, 
as shall be explained. 

First, a prestored trajectory cannot optimize fuel con- 
sumption or a similar performance measure under actual 
operating conditions. Optimum trajectories depend 
significantly on aircraft gross weight, wind and temperature 
profiles, and the initial state of the aircraft. These variables 
cannot be predicted with the required precision prior to 
takeoff. To prestore optimum trajectories for each of the 
conditions likely to be encountered would result in an im- 
possible large memory requirement. Therefore, prestored 
trajectories must necessarily represent a compromise in 
performance. 

Second, in existing systems, the pilot must fly the aircraft 
manually from its current position to the starting point of the 
trajectory. This flight segment is known as the capturing 
maneuver. Three-dimensional curved trajectories can be 
difficult to capture manually, and, if the trajectory also in- 
cludes a specification of landing time, as is the case in four- 
dimensional (4-D) guidance, the capturing maneuver cannot 
be done by the pilot without computer assistance. Therefore, 
the capturing maneuver, because of its variability, can only be 
generated by onboard trajectory synthesis. 
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Third, aircraft in high-density airspace are usually con- 
trolled by air traffic control vectors and during this period 
cannot follow a prestored flight path. Synthesis of a 
trajectory can only begin after the aircraft has received its 
final vector and has been cleared for approach. But the initial 
position of the aircraft at that time varies between ap- 
proaches; thus, trajectories require onboard synthesis. 

An initial design of a 4-D guidance system embodying the 
concept of onboard trajectory synthesis, including an ad- 
vanced capture law, was previously developed and flight 
tested onboard a conventional aircraft equipped with ad- 
vanced avionics.^ In the design described here, a new 
algorithm for generating horizontal capture trajectories has 
been implemented and vertical and speed profiles are syn- 
thesized using simplified aero/propulsion performance 
models of the aircraft. This method results in profiles that are 
more fuel efficient than those of earlier design. Design of the 
control law for tracking the synthesized trajectory is based on 
the linearized perturbation guidance approach. Since the 
perturbation equations are aircraft configuration dependent, 
gain scheduling is used in the feedback law. 

The Augmentor Wing Jet STOL Research Aircraft 
(AWJSRA) was chosen as the test vehicle for this concept. 
This type of powered-lift aircraft is highly cost-sensitive to 
operational procedures in the terminal area. It also exem- 
plifies particularly well the unique problems of powered-lift 
aircraft, namely, high fuel consumption in the STOL mode; 
dependence of both lift and drag on thrust; and an excess of 
controls over the minimum number needed to determine path 
and speed. These factors suggest that trajectory optimization 
could greatly increase the operational efficiency of the air- 
craft. Implementation of this concept was facilitated by the 
existing installation of an advanced avionics system onboard 
the aircraft. 


Equation (4) specifies the energy rate as a function of the 
difference between thrust and drag, subject to the constraint 
that lift equals weight. Thrust and drag are, in turn, functions 
of the controls producing forces in the flight -path direction; 
namely, throttle ir, flap angle nozzle angle (vectored 
thrust), and angle of attack a. Equation (5) determines the 
relationship between flight-path angle and deceleration for 
the energy rate calculated from Eq. (4). Equation (5) indicates 
that, in particular, a given energy rate may be utilized to fly at 
flight-path angle y„ with constant airspeed, or to fly at zero 
flight-path angle with acceleration dK^/dL An infinity of 
other combinations of and dK^/d/ can also be chosen to 
yield the same energy rate. This makes possible a simplifying 
dichotomy in the trajectory synthesis; namely, at any time the 
desired energy rate is selected first by choice of appropriate 
controls and then the linearly related quantities of y^ and 
dK^/d/ are selected to generate the specifics of the flight path. 

Since the STOL aircraft studied in this paper has four 
controls to achieve a specified energy rate and to maintain lift 
equal to weight, there is an excess of two controls over the 
minimum number needed for a simultaneous solution to Eq. 
(4) and the constraint L=fV, These two extra degrees of 
freedom in the controls are exploited to minimize power 
setting and, therefore, fuel flow at every energy rate. This 
optimization problem is restated in equivalent form as the 
maximization of energy rate for a given power setting: 

= max(T-D)/fV (6) 

r.Qt.iy 

Constraint: L(Tr,y,a,Sy) = IV 

The maximization must obey various inequality constraints 
on the controls: 


Energy Rate Model and Selection 
of Reference Controls 

An energy rate model of aircraft performance has been 
found to yield a compact and sufficiently accurate 
representation of performance for terminal area trajectory 
synthesis. In this section, a performance model based on 
energy rate is derived and then applied to determine the 
optimum reference controls for synthesizing trajectories. 
Consider the standard expression for energy rate written as 

d£ _ (r-D)K, 
d/ W 

where 

E^h^U/2g)Vi 

with constraint L=W (Ref. 3). It is assumed throughout this 
paper that flight-path angles are small such that cosy^ * 1 and 
sin » 7 a . Furthermore, it is assumed that flight-path angle 
rates are so small that their effect on lift is negligible. Dif- 
ferentiation of Eq. (2) with respect to time gives an equivalent 
expression for energy rate: 


( 1 ) 

( 2 ) 


d/ d/ g ^ dr 


( 3 ) 


Equations (1) and (3) can be nondimensionalized by dividing 
them both by Vg, The resulting quantity on the left side, 
(1/K^) (dE/dt ) , is defined as the normalized energy rate E„, 
or energy rate for short. By using the relation (dh/dt) * Valtat 
the two relations for E„ become 


En 


= (T-D)/W 
g dt 


( 4 ) 

( 5 ) 


- 10.5 deg^aS 19.5 deg 
6deg^»^s 100 deg 

5 . 6 deg ^ 6y ^ ( k;, ) [ flap placard ] 

In addition, a lift or maneuver margin must be satisfied at 
every point to guarantee sufficient normal force for changing 
the flight path. Pilots familiar with this aircraft specify that at 
least 0.4 g of normal acceleration must be attainable at any 
time by an increase in the angle of attack alone. 

The use of Eq. (6) results in the selection of the controls that 
yield the maximum attainable energy rate at each thrust 



Fig. 1 Energy rale diagram for STOL aircraft: 38,000 lb, sea 

level 59 *F. 


with constraint/. = W. 
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setting. This insures the efficient use of thrust at any energy 
rate that requires more than the minimum thrust. But energy 
rates more negative than those attainable by Eq. (6) are also 
of interest. Such negative energy rates must occur at the 
greater of the minimum or idle thrusts required by the 
maneuver margin. At a particular airspeed, a decrease in the 
energy rate below the minimum attained through Eq. (6) can 
be effected by increasing the vectored thrust angle v and/or 
the flap angle 5^, The third control, angle of attack a, is 
needed to satisfy the constraint L = 1V. The two degrees of 
freedom in the controls can be exploited to minimize noise 
exposure along the ground track. Noise under the aircraft is 
known to increase as the nozzles producing the vectored thrust 
are turned downward. Therefore, a further decrease in energy 
rate is achieved by first increasing flap angle until it reaches its 
limit or placard value and only then by increasing nozzle 
angle. 

The result of applying these procedures to the AWJSRA is 
shown in Fig. 1 for a weight of 38,000 lb, sea-level altitude, 
and standard temperature. The figure gives the envelope of 
energy rate vs indicated airspeed with throttle, flaps, and 
vectoring nozzle as parameters. Angle of attack is not plotted 
to avoid cluttering the figure. At any airspeed, the £„max 

curves define the range of permissible energy rates. The 
optimum controls for a given airspeed and energy rate are 
determined by interpolation between contours of constant 
controls. For example, at an airspeed of 105 knots and 
£„ = -0.17, the optimum controls are found to be 5^^ = 26 
deg, y = 6 deg, and ir = 84% (point A, Fig. 1). Angle of attack 
(not shown) is 8.4 deg. Maximum energy rate with minimum 
thrust occurs at 112 knots (point B) and corresponds ap- 
proximately to (L/D)^„ = 10. 

It should be noted that the force-producing controls in this 
experimental STOL aircraft have unusual characteristics that 
account for the relative complexity of Fig. 1. Throttle affects 
both lift and drag at all speeds, but the effect on lift is greatest 
in the STOL regime below about 80 knots. The thrust 
magnitude produced by the vectoring nozzle, referred to as 
the hot thrust, is also controlled by the throttle and accounts 
for about 60*7o of the total thrust produced by the two 
engines. The remaining 40*76 of the thrust, which is the cold 
thrust produced by the fans, energizes the augmentor wing to 
increase lift at STOL speeds. 

The relationship between the controls and the energy rate is 
revealed more clearly in Fig. 2 at the example airspeed of 105 
knots. Many such plots at various airspeeds would be required 
to illustrate the complete dependence of the controls on 
energy rate. As the energy rate decreases below its maximum 
value of 0.28, throttle decreases nearly linearly until idle 
throttle is reached. In this interval, flaps increase only 
slightly, while nozzle angle remains at minimum and angle of 
attack increases. At more negative energy rates, flaps become 
the dominant control until they reach the placard value of 40 
deg at this airspeed. Angle of attack decreases sharply as flap 
angle increases. Finally, nozzle angle increases toward its 



Fig. 2 Optimum controls as function of energy rate at 105 knots; 
38,000 lb, sea level 59*F. 


maximum value of 100 deg as the energy rate decreases 
toward its negative limit of -0.3. 

In the flight implementation of the algorithm, four 
diagrams, as shown in Fig. 1, are utilized: two for sea-level 
altitude at weights of 38,000 and 48,000 lb, and two others for 
5000-ft altitude at similar weights. Experience indicates that 
these provide sufficient data to adequately interpolate the 
controls. Each diagram requires 124 words of memory in the 
airborne computer. The small circles in Fig. 1 indicate the 
locations of points that are stored. The energy rate data are 
also corrected for deviations from the standard temperature 
profile. Correction is done by computing a thrust setting 
corrected for temperature deviations. 

Synthesis of Complete Profiles 

In the preceding section the criteria of fuel conservation 
and noise reduction were used to determine the four reference 
controls of throttle, nozzle angle, flap angle, and angle of 
attack as functions of the energy rate. This approach replaced 
the problem of selecting four control variables with the 
simpler problem of selecting a single, equivalent variable, 
namely, the energy rate. In this section, we make use of the 
energy-rate variable in generating efficient terminal area 
trajectories. 

The problem of terminal area trajectory synthesis can be 
stated as the specification of rules for flying an aircraft with 
initial state vector [Xf, Y,, //,, to a final state vector 

[Xp Yp hp Hp . To be of practical interest, such rules 
must generate efficient and flyable trajectories connecting 
various initial and final state vectors. By specifying a per- 
formance criterion such as fuel consumption, we can fit this 
problem into the framework of optimal control theory. 
However, the difficulty of solving an optimal control problem 
characterized by a five-element state vector makes this ap- 
proach computationally impractical for in-flight im- 
plementation. Following Ref. 4, we have adopted the sim- 
plifying procedure of separating the synthesis problem into 
two essentially independent problems. 

The first problem consists of synthesizing the horizontal or 
2-D trajectory. References 4 and 5 give algorithms for 
computing near -minimum-distance 2-D trajectories as a 
sequence of an initial constant radius turn, straight flight, and 
a final constant radius turn, where the turn radii are chosen so 
as to avoid exceeding a specified maximum bank angle at the 
maximum groundspeed encountered in each turn. The 
algorithm used in this onboard computer implementation is 
based on a simplified derivation which resulted in a 
significant reduction in computing time compared with the 
methods given in the references. The derivation can be found 
in the Appendix. Figure 3 illustrates two of the four types of 
horizontal trajectories that can occur; the two not shown 
differ only in the orientation of the turns. The algorithm 
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Fig. 3 Examples of minimum distance, constant turn radius, 
horizontal capture trajectories to a capture point Pt on final ap- 
proach: a) turns in same direction; b) turns in opposite direction. 
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computes all feasible trajectories, of which there are at least 
two and at most four, and chooses the one with the shortest 
path length. Note that the terminal point lies on an extension 
of the runway centerline and that the final heading angle is 
equal to the runway heading. The final point should be chosen 
as close as possible to the touchdown point consistent with 
safe operational practice. For a STOL aircraft, the minimum 
distance is about 1 n.mi. 

The second problem, solved after the horizontal trajectory 
has been computed, consists of synthesizing efficient speed 
and altitude profiles which match the initial and final speeds 
and altitudes and respectively. 

The horizontal distance of the trajectory df^, a known 
quantity computed in the previous step, adds a third boundary 
condition to be satisfied by the profiles. While this three-state 
optimal control problem is much simpler to solve than the 
original five-state problem, it is still too complex for onboard- 
computer implementation. A simpler algorithm was, 
therefore, developed that generates near-optimum speed- 
altitude profiles by matching the general characteristics of 
optimum fuel and noise trajectories studied in Refs. 6 and 7, 
respectively. We briefly explain the rationale for this 
algorithm with reference to descent, which is the most dif- 
ficult case. 

It was found in Ref. 6 that the descent portion of a 
minimum-fuel descent trajectory is characterized by a delay in 
the start of the energy decrease as long as possible, consistent 
with meeting end constraints of speed and altitude. Fur- 
thermore, the energy change consists initially of descent to the 
final altitude at near-constant indicated airspeed followed by 
deceleration in level flight. Most of the energy change takes 
place at minimum throttle, as one might expect for minimum 
fuel flight. Minimum-noise descent profiles computed in Ref. 
7 are similar in that they also delay the start of energy 
decrease as long as possible, but approach the final altitude in 
a steep descent to maximize the aircraft’s altitude above the 
ground near the runway. This means that the deceleration to 
the final airspeed takes place before the start of descent or 
during the early portion of the descent. Thus, the two types of 
descent profiles differ primarily in the way they proportion 
the use of available energy rate to decrease altitude and 
airspeed. 

To facilitate the synthesis of such profiles, a family of 
decreasing (and by extension, increasing) energy profiles, 
which include the two types described as special cases, is 
defined by two parameters, o and e. The first parameter, a, 
selects the fraction of minimum /maximum available energy 
rate, (^nmax) be used for decreasing/increasing 

energy. The values of £„min ^ud ean be read from Fig. 1 
at each indicated airspeed. The second parameter, c, deter- 
mines the fraction of the selected energy rate to be used for 
deceleration/acceleration. Then, for particular choices of a 
and e, the energy rate, airspeed, flight-path angle, altitude, 
and horizontal distance are computed as follows: 



(7) 

V,=gtE„ (0rSeS7) 

(8) 

7, = (/-0£„ 

(9) 

ii=y.yo 

(10) 

S= K.COS7, + 

(11) 


where is the along-track component of windspeed. Note 
that Eqs. (7-9) are consistent with Eqs. (4) and (5) for all 
values of a and c. Decreasing/increasing energy profiles are 
generated by integrating Eqs. (8), (10), and (1 1) for particular 
choices of a and e. 

To illustrate the effect of the parameter e on the 
descent/deceleration profiles, assume £„=-0.13, in- 




Fig. 4 Effect of € on speed and altitude profiles, with E„ = - 0.13. 


dependent of speed, and let the airspeed to be achieved at 
touchdown be 100 ft/s. To achieve the desired boundary 
conditions, Eqs. (8), (10), and (11) are integrated in backward 
time starting with the speed and altitude at touchdown. The 
resulting airspeed and altitude profiles are plotted as a func- 
tion of distance to touchdown in Fig. 4 for € = 1, 0.5, 0.0. The 
profile for € = 1 is seen to approximate the minimum fuel, for 
c = 0, the minimum noise descent and for c = 0.5, a com- 
promise between fuel and noise minimization. 

To minimize fuel consumption or noise, changes in energy 
should be made at maximum rate when the aircraft enters the 
powered-lift region of 90 knots and below. This is ac- 
complished by setting a to unity and thereby following the 
£„nijn contour during descent and deceleration. However, for 
the aircraft under study this can yield energy rates too 
negative for safe operation. A limit less than one is also 
necessary to reserve energy rate for perturbation control. A 
practical upper limit on a is about 0.9 for the AWJSRA. In 
the flight implementation, the two profile parameters are 
keyboard entries that allow the pilot to choose values ap- 
propriate for each landing approach. In addition, the pilot 
can specify the maximum deceleration and descent angles via 
keyboard entry. The maximum safe deceleration for this 
aircraft is limited to about 0.06 g by the maximum rate of 
which flaps can be extended. The synthesis algorithm is 
configured to decrease a below its limit if that is necessary to 
satisfy these constraints. 

The backward time integration described above generates 
an increasing (in backward time) energy profile starting at the 
desired final speed and altitude. To complete the synthesis of 
the descent trajectory, we still need rules for matching this 
profile to the initial speed and altitude of the aircraft. The 
freedom of the aircraft to maneuver in altitude is restricted by 
air traffic control as well as passenger comfort con- 
siderations. Thus, as an aircraft approaches a terminal area, it 
is generally not allowed to climb above its initial approach 
altitude for the purpose of optimizing the approach trajec- 
tory. The aircraft must hold this altitude until starting the 
final descent. However, while flying at altitude /?,, it may 
change to a new airspeed, called the terminal area sp>eed, 
which can be higher or lower than the initial speed . Unless 
specified by the pilot via keyboard entry, it is chosen to 
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minimize fuel use per unit distance, and is 140 knots for this 
aircraft (it would be 220-250 knots for conventional jet 
transports). 

The various rules contained in the preceding two 
paragraphs can now be combined to yield the complete 
algorithm. The synthesis begins with the backward time in- 
tegration from final conditions using the specified a 

and e. If the altitude reaches its target value of /», before the 
airspeed reaches its target value of we set e = 1 and then 
continue the backward time integration until the airspeed has 
also achieved its target value. When setting €= 1, the flight- 
path angle is forced to zero and the energy rate is used entirely 
for accelerating (in backward time) toward On the other 
hand, if the airspeed reaches its target value before the 
altitude does, we set e = 0. This stops the airspeed change and 
uses the energy rate entirely for increasing the altitude toward 
its target value of h^. When the second and last variable 
reaches its target value, we set (j = 0, i.e., E„ =0, thus com- 
pleting the backward time integration. Next, we begin a 
forward time integration to get the distance required to 
change speed from to with c = 1 . Let the distances for 
the backward and forward integrations be df, and d^, 
respectively. A valid trajectory has been generated if the 
cruise distance d^., computed from 

dc = d,-d,-dj ( 12 ) 

is nonnegative, i.e., If d^ is negative, the synthesis has 

failed because the aircraft is too close to the capture point Pj. 

Figure 5 illustrates the various segments of an approach 
trajectory synthesized by the algorithm. As before, we assume 
for simplicity that = — 0. 13, a constant. Other parameters 
defining the problem are indicated in the figure. Note that the 
initial descent at = -7.5 flattens to 7^ = - 3.75 to allow the 
aircraft to decelerate. The reference controls for this 
trajectory can be interpolated from Fig. 1 . 

The airspeed deceleration is corrected for known wind 
shears, which are computed from a knowledge of V^{h), if 
available. The wind shear correction factor is 

AK,= -(dK,/d/i)K,7. 

and is added to the right side of Eq. (8) to obtain ihe corrected 
airspeed rate. Furthermore, the reference controls are 



Fig. 5 Example of synthesized STOL approach trajectory . 


corrected for the effect of the bank angle used in flying a turn 
by interpolating the controls with lift equal to an aircraft 
weight multiplied by the load factor \/cos<t>. Integration step 
size varies during synthesis. During decelerations or ac- 
celerations it is 1 s, while during altitude changes at fixed 
speed it is 5 s. Total time for synthesizing a complete 
trajectory consisting of a horizontal trajectory similar to the 
ones shown in Fig. 3 and a speed /altitude profile similar to 
the one in Fig. 5 is about 2 s on the particular airborne 
computer used in the flight tests. When the trajectory syn- 
thesis is time-shared with navigation and other necessary 
computations, the computing time increases to about 6 s. 

Perturbation Guidance Law 

Perturbations of the aircraft states from the reference states 
are used in the guidance law to generate perturbation controls 
which are added to the reference controls in order to null 
errors in airspeed, altitude, and crosstrack position. The 
feedback states in the guidance law also include crosstrack 
error rate and flight-path angle, as well as the integrals of 
airspeed and altitude errors. The latter two are used to reduce 
speed and altitude bias errors caused by inaccuracies in the 
stored energy rate data and errors in the estimates of wind and 
temperature profiles. 

The controls are throttle, nozzle, pitch, and roll angles. 
Flaps are not used as perturbation controls because of their 
relatively low rate limit and an operational constraint that 
flap motion be monotonic during an approach. The flap 
command is simply the reference value at each ground track 
position limited to the placard value at the current airspeed. 

Lateral perturbation control is essentially uncoupled from 
the longitudinal mode and is accomplished through a roll- 
angle command to the roll-command autopilot. This com- 
mand is of the form 

<t>^ = <t>r^k^yAy-^k^^Ay 

where <t>r is the reference roll angle, and Ay and Ay are the 
crosstrack error and error rate, respectively. The two gains 
were chosen to provide a well-damped response and control 
activity compatible with the noise characteristics of the 
navigation system. 

Longitudinal perturbation control for correcting airspeed 
and altitude errors is difficult because the reference controls 
generated by the energy rate schedule of Fig. 1 often lie on a 
constraint boundary and therefore cannot be perturbed freely 
in both directions. The two controls that are often constraint- 
limited during a fuel-conservative approach are throttle x, 



Fig. 6 Constraint boundaries for perturbation controls. 
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and nozzle angle v. Some insight into this problem can be 
obtained using data from the energy rate schedules. Figure 6 
shows the energy rate envelope from Fig. 1 with the minimum 
reference nozzle and throttle constraint boundaries. These 
boundaries divide the envelope into four regions: I, where v 
cannot be reduced; II, where neither x nor v can be reduced; 
III, where x cannot be reduced; and IV, where x and p are free 
to move in either direction. The combinations of controls 
available for increasing and decreasing E„ in each region are 
indicated in the figure. Note that in region I, the nozzle could 
be used as a additional control variable for decreasing energy 
rate. However, this variable is not used because throttle and 
pitch provide adequate control of flight-path errors in this 
region. In region IV, the minimum reference throttle is above 
idle and is determined by the maneuver margin constraint. At 
each airspeed in this region the negative throttle perturbation 
that can be added to the reference throttle to yield the com- 
manded throttle is limited to -2% for safety reasons. 
Positive and negative throttle perturbations are further 
limited so that the commanded throttle x^ falls into the engine 
operating range, 84 Vo ^ x^ 96 Vo. 

The perturbation equations and the perturbation control 
law can be written in state vector notation as 

dx/dt=Fx-\-Gu u = Kx 

where 

x= (AF,A7,AA,jAFdr,iAAdO ^ 

«=(Ax,A^,Ai/)^ 

The delta quantities are the perturbations from reference 
values, i.e., AV= - etc., where is the aircraft and 
Vgr the reference true airspeed, respectively. The commanded 
controls are the sum of reference and perturbation controls: 

W(.= (x^ + Ax,0^ + + A*^) 

For a powered-lift STOL aircraft, such as the one used for 
these flight tests, the values of Fand G are strongly dependent 
upon airspeed and energy rate and are thus time-varying along 
a trajectory. Quadratic optimal synthesis* would therefore 
yield time-varying gain matrices that are also functions of the 
reference trajectory. But it is neither practical nor necessary to 
implement a complex, reference-trajectory-dependent gain 
matrix in order to achieve adequate control system per- 
formance in this case. 

The design procedure employed here began by first com- 
puting optimum gain matrices at various operating points in 
the control region diagram (Fig. 1) using fixed values of Fand 
G. The analysis of these gain matrices showed the strongest 
dependence on airspeed, reference nozzle angle, and reference 
flaps. Sensitivity of the closed-loop eigenvalues to changes in 
several of the gains was low, allowing those to be set to zero or 
held constant throughout the operating region. It was possible 
to fit the variable gains with relatively simple functions of 
reference airspeed, nozzle angle, and flap angle. This method 
resulted in the following gain matrix: 



Fig. 7 Simplified flow diagram of guidance system. 


where is in units of feet per second. Extensive computer 
calculations have verified that the closed-loop eigenvalues of 
this system have damping factors of 0.707 or greater and real 
parts less than - 0.05/s at all operating points. These 
characteristics provide adequate tracking performance. When 
operating in region I of Fig. 6, the last row of /T is set to zero 
since nozzle angle is not used for control. In regions II and 
III, throttle perturbations are limited to positive values, while 
in region II, nozzle perturbations are limited to positive 
values. In region IV, each control moves freely but negative 
throttle perturbations are limited to -2% rpm, as previously 
explained. Control limiting can reduce the effectiveness of 
integral feedback of speed and altitude. Some design con- 
siderations for these integral feedback loops are given in 
Ref. 9. 

The throttle and nozzle angle perturbations generated by 
the control law will generally be of opposite sign, because the 
elements of the first row of K all have opposite sign of the 
third-row elements. Thus, even in region II, where throttle 
and nozzle perturbations are each limited to move only in the 
positive direction, they are not generally limited 
simultaneously. This implies that two controls, either throttle 
and pitch or nozzle and pitch, are free to move. Transient 
response studies using a nonlinear simulation of the aircraft 
and guidance system have shown that the control power is 
adequate to provide rapid and well-damped airspeed and 
altitude error responses in region II . 

Guidance Algorithm Overview 

A flow chart illustrating the integration of major functions 
within the guidance system is shown in Fig. 7. The pilot enters 
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into the guidance computer the coordinates, altitude, and 
landing airspeed to be achieved at the desired final position on 
the approach path (see Fig. 3) or, alternatively, he selects a 
capture waypoint with prestored coordinates. Trajectory 
synthesis can begin after the navigation system has computed 
the current position and velocity components of the aircraft. 
The first step in the synthesis process involves computing the 
horizontal trajectory parameters using the technique 
described in the Appendix. This step is always successful. In 
the second step, the altitude and speed profile are synthesized 
using the energy rate diagram of Fig. 1 in conjunction with the 
logic described in the section on Synthesis of Complete 
Profiles. This step is not always successful. For example, if 
the horizontal path computed in the first step is very short and 
the differences in speed and altitude between initial and final 
aircraft positions are large, a flyable trajectory along that 
path may not exist. During synthesis, such a failure is detected 
as a negative cruise distance in Eq. (12). A failure to 
synthesize is an unlikely event in landing approaches initiated 
several miles from the final point, the usual situation; if the 
failure occurs, however, the synthesis is repeated using up- 
dated position and velocity vectors, which aircraft motion has 
changed during the time the synthesis was in progress. The 
pilot can also fly the aircraft manually to a more favorable 
location (path stretching) for a successful synthesis. 

After a trajectory has been successfully synthesized, its 
parameters are stored and the horizontal path is displayed to 
the pilot on a map-like cathode ray display (see Ref. 2 for 
details of this device). The appearance of the trajectory on the 
display is also a cue to the pilot that a valid trajectory has been 
successfully synthesized. The stored and displayed trajectory 
is refreshed by repeating the synthesis every few seconds until 
the pilot engages the auto track mode. At that time, the last 
synthesized trajectory is frozen and its reference states and 
controls are regenerated in real time. These reference values 
are fed forward to the perturbation feedback law, which 
causes the aircraft to track the synthesized trajectory. Three- 
dimensional navigation data for computing errors between 
the aircraft and the reference trajectory are obtained from 
TACAN (TACtical Air Navigation System) or MODILS 
(MODular Instrument Landing System, an experimental 
Microwave Landing System), with automatic switching to the 
most accurate signal. 

It is important to note that once automatic tracking of a 
trajectory has begun, the stored trajectory is not refreshed, 
though this may be desirable if unmodeled wind or transients 
in navigation introduce large tracking errors. The real-time 
operating system has been configured to add this capability in 
future flight experiments. Also, a technique is incorporated in 
the flight software that compensates the trajectory for the 
changes in aircraft position occurring during the time (up to 6 
s) the trajectory is synthesized. 

Flight Test Results 

Figure 8 shows the major portion of various time histories 
for a straight-in flight test approach starting 7 n.mi. from 
touchdown, at 3000 ft altitude and 140 knots. The 
deceleration at 0.03 g begins in level flight at point A, 4.7 
n.mi. from touchdown. The descent begins at point B, 3.7 
n.mi. from touchdown. The entire approach trajectory was 
flown using TACAN for navigation. The TACAN station was 
located at the airport a few hundred feet from the runway 
centerline. Because of its favorable location, the TACAN 
station provided sufficient navigation accuracy for flying the 
approach automatically to within a half-mile of touchdown 
without switching to the higher precision MODILS as would 
normally be required. There was light-to-moderate turbulence 
and an average headwind of about 1 5 knots below 4000 ft as 
measured by a radar tracked weather balloon just prior to 
takeoff. However, the wind profile was not entered into the 
synthesis logic and thus constituted an unmodeled wind. 
Altitude errors, except near the pitchdown point, did not 
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Fig. 8 Flight test results, straight-ln approach, runway at 140 ft 
above sea level. 

exceed 35 ft and decreased to about 15 ft near the end. Speed 
errors during deceleration were less than 10 ft/s, and 
decreased to about 1 ft/s at the end. If allowance is made for 
the presence of turbulence, winds, and navigation system 
noise during the flight, these errors agree reasonably well 
with simulation results and are acceptably low. The control 
perturbation biases, evidently caused by modeling errors and 
the unmodeled wind, are larger than those seen in simulation, 
though they are not excessive. Nozzle bias during the middle 
of the deceleration averages about 25 deg. While this seems 
large, it should be noted that during this interval the throttle is 
at flight idle, where the effect of nozzle on energy rate is a 
minimum. On the whole, the control biases represent fairly 
small errors in the energy rate model. The flight test results 
can, of course, be used to improve the accuracy of the energy 
rate model of the aircraft. 

The crosstrack error at the end of deceleration, point C, 0.9 
n.mi. from touchdown, was measured by precision radar as 
80 ft. This error is an important criterion for determining how 
close point C can be placed to touchdown when navigating 
with TACAN. Pilots judged the action of the automatic 
control law as smooth and the trajectory synthesis technique 
as a convenient and effective tool for optimizing approach 
trajectories. 

The fuel consumption of this automatically flown 
trajectory was compared with that of a trajectory flown by a 
test pilot under simulated instrument fligh rule conditions. In 
order to provide a basis for comparison, the manually flown 
trajectory began from the same initial distance-to-touchdown, 
airspeed, and altitude as the automatically flown trajectory. 
The approach was made with the aid of a flight director 
system which displayed to the pilot lateral and longitudinal 
deviations from a straight -in 7.5 deg approach path. The fuel 
used for the automatic approach was 381 lb, while that for the 
manually flown one was 500 lb. Further simulations and 
flight tests are in progress to compare the fuel consumption 
for various approach trajectories, flight director designs, and 
wind conditions. 

Conclusions 

The automatic guidance system described in this paper 
achieves the dual goal of fully automatic flight and near- 
optimal fuel conservation through the technique of fast-time 
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onboard trajectory synthesis. This technique overcomes the 
performance limitations inherent in a stored precalculated 
trajectory by adapting the trajectory to the unique conditions 
encountered in each landing approach. The ability to adapt is 
crucial in the terminal area since the initial conditions for 
starting the approach and the wind and temperature profiles 
are not predictable with sufficient accuracy prior to takeoff. 
Synthesized profiles always delay the start of the descent and 
deceleration points as much as possible. Flight evaluation 
using a powered-lift STOL aircraft showed that an 
automatically synthesized approach saves approximately 120 
lb of fuel during the last 7 miles of the approach relative to 
one flown manually with only conventional flight director 
guidance. The design procedure described herein for a STOL 
aircraft is applicable with lesser computer complexity to 
conventional aircraft. The algorithm is also suitable for in- 
corporating in advanced flight management systems currently 
under development by industry. 


and therefore, since D and are perpendicular, 


Dm \D\ (A2) 

where, by definition, 

G- IQI =y/{XC,-XCj)^ + (YC2-YCf)^ (A3) 

It can be seen from Eq. (A2) that no real solution exists if 
G< \R2S2-R jSj I. When the turns are in opposite direc- 
tions, S/ = — S2 , and there is no real solution for 
G< (/? / +^2) » i.e*, if the circles intersect. On the other hand, 
for rotations in the same directions, 8^=82 and a real 
solution exists unless Q< 1/?^ -/?; 1 , i.e., unless one circle lies 
entirely within the other. From geometric construction it can 
be shown that there always exist at least two real solutions. 
From the definition of the radius vectors, one can write for 
the real solutions: 


Appendix 

This appendix derives the expressions for synthesizing 
horizontal capture trajectories for flying an aircraft from a 
given initial position and heading to a specified final position 
and heading. Figure 3 is used to explain the problem and 
define the variables. The turns are arcs of the circles shown in 
the figure and the straight portion of the trajectory must be a 
line tangent to both circles. Since the initial and final turns 
may be either clockwise or counterclockwise, there are four 
possible combinations of turning directions — two with the 
initial and final turns in the same direction, and two in op- 
posite directions. Figure 3 illustrates one solution of each 
type. If a given pair of circles is entirely separate, i.e., no part 
of one circle lies within the other, it is possible to draw four 
tangent lines between the pair. However, vector D along the 
tangent line from the initial to the final circle coincides with 
the direction of rotation at both tangent points for only one of 
the four tangent lines as shown in the figure. 

In the figure, the final position and the origin of the 
coordinate system are located on the runway centerline. 
However, the derivation is for arbitrary locations. Fur- 
thermore, all variables are defined so that the derivation 
applies to all possible combinations of turning directions. 

Figure 3a is for the case where both turns are in the same 
direction and the tangent vector D does not cross G* while in 
Fig. 3b the turns are in opposite directions and D crosses Q. 
Initially, the aircraft is at {X^, T,) in some inertial Cartesian 
coordinate system with heading //, defined as positive 
clockwise from the X axis, and Vj is a unit vector in the 
direction of the velocity. The vector distance from (X^, Y^) to 
the center of the turn is given by UjRj, where Rj is the radius 
of turn and My a unit vector normal to Vj and positive to the 
right of Uy. Therefore, the vector from (A',, T,) to the center 
(A'Cy,yCy) is RjUj for a right turn and —RjUj for a left 
turn. The directions of turn are accounted for by writing the 
radius vector as /?ySyi/y, where Sy = + 1.0 for right turns and 
S; = — 1.0 for left turns. Similarly, the direction of the final 
turn is denoted by S^. 

The aircraft moves along the circle from (Ay, Ty) to the 
tangent point ( ) which has a radius vector ^ySy m^. The 
tangent vector from (A^, initial turn to 

(Aj, y^) at the beginning of the final turn is D, The radius 
vector at (Aj,Tj) is but since m^ and lij must be 

normal to D, m^ = lij . Likewise, the headings H2 and at 
the two tangent points are equal. The final turn ends at 
(Xp Yj) with heading Hf and radius vector 

Using this notation, we can write 


and 


R,u,S,= 


/-/?ySySin//y\ 
\ /?ySyCOS//y / 


R,u,S, 



(A4) 


(A5) 


Equating Eqs. (A4) and (A5) gives 

ACy = Ay -/?ySySin//y YC j = Yj +/?ySyCOS//y (A6) 

Similarly, 

AC^ =A^-/?jS^sin//^ = Tyr-Z^^S^cos//^ (A7) 

The radius vectors at the tangent points can be used in the 
same manner to compute the components of A^ and Aj 


X2 = ACy + /? / Sy sin//^ (A8a) 

Y2 = yCy -Rj 8 ^qo%H 2 (A8b) 

A3 = AC3 + /?3S3Sin//3 (A9a) 

= YC2-R282C0%H2 (A9b) 


Subtracting Eq. (A 8a) from (A9a) and Eq. (A8b) from (A9b) 
gives the components of D\ 

A3 - A3 = AC3 - ACy + (R282-RJ8J ) sinH2 

(AlO) 

Tj - y^ = YC2 - YCj - {R282-Rj8j)cosH2 

Another expression for the components of D is: 

A3-A3 =Dcos //3 Y2-Y2=Ds\nH2 (All) 

Equating the corresponding pairs in Eqs. (AlO) and (All) 
gives 

DCOSH2 = (AC3 - ACy ) + {R282 -/?ySy )sin//3 

(A12) 

Ds\nH2 = ( yc^ - yCy ) - (/?3S3 -/?ySy ) cosHj 


D-\-R2U2^2^^1 ^2^1 Q 


Equations (A 12) can be solved for the tangent of H2 : 


or 


tan//^ 


(yC3-yCy)D-(/?353-/?ySy)(AC3-ACy) 

(AC3 - ACy )D+ (R282 -Rj8j) ( YC2 - yCy ) 

(A 13) 


Q— D U2 (R282 — /?ySy ) 


(Al) 
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Equations (A6-A9) and (A13) completely specify a capture 
trajectory for any combination of Sy and S^. However, the 
length of the trajectory is also needed in order to determine 
which of the feasible trajectories gives the minimum distance. 
The first turn angle is 

TRj = (Hj-H,)+27tC^Sj 

where 

CO if 

^ 1/ if Sj{H2-Hi)<0 (A14) 

and the second turn angle is 

TR2={Hf-H2)^2TrC2S 

where 

CO if S2{H^-H2)^0 
^ 1; if S2{H^-H2)<0 (A15) 

Finally, the total length of the capture path is 

d^= (A16) 

The algorithm computes the length of trajectories for all 

feasible pairs of S; and and then picks the shortest length 
trajectory. A FORTRAN listing and additional details of this 
algorithm, including captures with three turns, can be found 
in Ref. 10. 
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